Empirical Fit to Precision Inclusive Electron-Proton Cross Sections in the Resonance 

Region 

M.E. Christ£] 

Hampton University, Hampton, Virginia 23668 

P.E. BostecQ 

Jefferson Lab, Newport News, Virginia 23606 
(Dated: May 6, 2010) 

An empirical fit is described to measurements of inclusive inelastic electron-proton cross sections 
in the kinematic range of four-momentum transfer < Q 2 < 8 GeV 2 and final state invariant 
mass 1.1 < W < 3.1 GeV. The fit is constrained by the recent high precision longitudinal and 
transverse (L/T) separated cross section measurements from Jefferson Lab Hall C, un-separated 
Hall C measurements up to Q 2 ~ 7.5 GeV 2 , and photoproduction data at Q 2 = 0. Compared to 
previous fits, the present fit covers a wider kinematic range, fits both transverse and longitudinal 
cross sections, and features smooth transitions to the photoproduction data at Q 2 = and DIS data 
J>^' at high Q 2 and W. 



I. INTRODUCTION 

Knowledge of the inclusive electron-proton cross section in the nucleon resonance region is important input for 
many research activities in nuclear and particle physics. Classic examples include calculations of radiative corrections 
to cross sections and extractions of spin structure functions from asymmetry measurements. In addition, the latter 
requires knowledge of the separated longitudinal and transverse photo-absorption cross sections. More recently, the 
electron-proton resonance region cross section has been used in constraining the vector coupling in models [H Q of 
low energy neutrino-nucleon cross sections. This is important since the quality of low energy neutrino-nucleon cross 
section models will become one of the largest uncertainties in the extraction of neutrino oscillation parameters from 
future long-baseline experiments. 

Having the best possible fit to the inclusive electron-proton cross section is important for all the endeavors described 
(■ — \ above. In this paper we will describe a fit to precision measurements of inclusive electron-proton cross sections in 
£f) • the resonance region for Q 2 < 7.5 GeV 2 . Among the advantages of this fit over previous resonance region proton 
fits 0, 0] is the constraint on the individual transverse and longitudinal cross sections afforded by the use of the recent 
I . high precision L/T separated cross section measurements [5j, |6[ from Jefferson Lab Hall C, a smooth transition to the 
j photoproduction point, and use of threshold-dependent Breit-Wigner forms for all resonances. It needs to be stressed 
that while the fit form utilized is, in general, physically motivated, the focus of the current work is to provide the 
best description of the inclusive proton data and not to determine the masses, widths, transition form factors, and 
branching ratios of the produced resonant states. These are best determined from exclusive meson production data. 
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II. INCLUSIVE ELECTRON SCATTERING FROM THE PROTON 

In terms of the incident electron energy, E, the scattered electron energy, E , and the scattering angle, 9, the 
absolute value of the exchanged 4-momentum squared in electron-proton scattering is given by 

Q 2 = (-q) 2 =4EE sin 20 -, (1) 

and the mass of the undetected hadronic system is 

W 2 = M 2 + 2M p v-Q 2 , (2) 
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with M p the proton mass and v = E — E' . In these expressions we have neglected the electron mass which is negligible 
for the kinematics studied. 

In the one-photon exchange approximation, the spin-averaged cross section for inclusive electron-proton scattering 
can be expressed in terms of the photon helicity coupling as 

^ = r [a T (W 2 , Q 2 ) + ea L (W 2 , Q 2 )] , (3) 
where a? is the cross section for photo-absorption of purely transverse (longitudinal) polarized photons, 



aE'(W 2 - M 2 ) 
(2n) 2 Q 2 M p E(l - e) 



r = v — (4) 



is the flux of virtual photons, and 



l + 2(l + ^)tan 2 - 



(5) 



is the relative flux of longitudinal virtual photons. Since T and e are purely kinematic factors, it is convenient to 
define the reduced cross section 

.^^ = .,(^g 2 ) + ^(^ 2 ,g 2 ). (6) 

All the hadronic structure information is, therefore, contained in ox and ox, which are only dependent on W 2 and 
Q 2 . 

In light of the current discrepency between the proton elastic form factors extracted from unpolarized scattering 
to those extracted from polarization observables [7j, it is worth commenting on the assumed validity of the one- 
photon exchange approximation. While, the most commonly accepted explaination for resolving this discrepency is 
the existance of two-photon (2-7)exchange contributions to the elastic scattering, the available data indicate that the 
reduced cross section is still quite linear in e over the ranges measured. Indeed, a recent search for non-linearities 
in the combined elastic and inelastic data which separated the longitudinal and transverse cross sections (including 
a significant portion of the resonance region and DIS data which is fit here) found no significant (i.e. less than 1% 
typically) evidence for non-linearities [g] . Hence, although it is impossible to rule out the existence of 2-7 contributions 
to the inelastic cross sections from this data, the existence of such contributions does not affect our ability to provide 
a reliable fit to the measured cross section data within the one-photon approximation. It is only the interpretation of 
the separated structure functions which could be incorrectly interpreted if significant 2-7 contributions exist. 



III. DESCRIPTION OF THE FIT 



In the past, fits to resonance region cross sections have typically not included data covering a large enough range in 
e and/or of high enough precision to constrain both ox and ox- Therefore, such fits have relied on an educated guess 
for the ratio R — o\l/<7t in order to extract data on ot from the measured reduced cross sections, and it was the 
extracted ot that was consequentially fit. This guess has typically been to use a fit to deep inelastic scattering (DIS) 
data on R extrapolated into the resonance region. The assumption here is that all the resonant structure cancels 
in R. However, it is now clear from the recent Jefferson Lab L/T data 5] that the resonant structure in o^ differs 
significantly from that in <jt, resulting in W and Q 2 — dependent variations of up to 0.1 in R. 

Utilizing the new precision data, the reduced cross section was fit by parameterizing ctt(W 2 ,Q 2 ) and <7l(W 2 , Q 2 ) 
and then minimizing the difference of ay constructed from Equation [5] with respect to the data. The cross section 
parameterization contained 75 free parameters which were determined from the fit: 7 for the resonance masses, 7 for 
the resonance widths, 25 to describe the Q 2 dependence of the transverse transistion form factors, 18 to describe the 
Q 2 dependence of the longitudinal transistion form factors, 10 to describe the non-resonant contribution to cft, 7 to 
describe the non-resonant contribution to ox, and a single damping parameter for the delta resonance. A detailed 
descripton of the parameterization is given later in this section. 

In order to avoid local minima and to speed up the convergence the starting parameters for ot and ox were 
determined by first fitting the approximate separated cross sections independently. To accomplish this the following 
procedure was utilized: first, the i?igg8 [9] fit to DIS data on R was used to extract ot and ox from oy; second, 
the extracted data on ot and ox were fit independently, and finally, R determined from the new fits was used to 
extract ot and ox for the next iteration. Several iterations were performed to ensure that reasonable convergence was 
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obtained. The fit forms for the resonant and non-resonant contributions to the cross sections were determined from 
this procedure for both ot and and were then used in a single fit to ay in which all transverse and longitudinal 
parameters were allowed to vary simultaneously. 

The fit form used to describe the separated resonance region photo- absorption cross sections was based on the 
following ansatz: 1) the cross section is the incoherent sum of contributions from resonance production (a R ) and a 
non-resonant background (a NR ), 2) the resonant cross section can be described by threshold-dependent relativistic 
Breit-Wigner forms with Q 2 -dependent amplitudes for each resonance, and 3) the non-resonant background varies 
smoothly with W 2 . Therefore, 



NR 
T,L> 



(7) 



with resonant contribution 



4.l(W 2 ,Q 2 ) = WJ2BW^ L (W 2 ) ■ L4^, L (Q 2 )] 2 



(8) 



The Breit-Wigner form utilized is 



BW l = 



KiKf 



1 i 1 i 



with 



KK cm [{W 2 - M 2 ) 2 + (MiTf 1 ) 2 } ' 
W 2 - Ml 

K -- 



2M„ 



(9) 



(10) 



K c 



W 2 - M 2 
ZW : 



(11) 



Ki = K 



M, 



and 



R cm = R c 



\Mi- 



(12) 



(13) 



Here, K and K cm represent the equivalent photon energies in the lab and center of mass (CM) frames, respectively, 
while Ki and Kf m represent the same quantities evaluated at the mass of the i th resonance, Mi. r*°* is the full decay 
width defined by 



-\tot 



(14) 



with f3j the branching fraction to the j th decay mode for the i th resonance and Tj the partial width for this decay 
mode. For single-meson decay modes the partial widths were determined from 



cm 

1 i 


21 + 1 


~(p™\ Mi ) 2 +X 2 ~ 


Pj 1 Mi 




(pf) 2 + x? 



(15) 



where Ti are the intrinsic widths of each resonance, the Pj m are meson momenta in the center of mass, I is the angular 
momentum of the resonance, and Xi is a damping parameter. In the present analysis it was found that a good fit 
to the data was obtained using Xi = 0.215 GeV for all the resonances except the P33(1232), where a smaller value 
of Xi = 0.1446 GeV resulted in a modest improvement. Overall the goodness of the fit was found to be relatively 
insensitive to the exact value of Xi. For the two-pion decay mode the partial widths were determined from 



Mi 



-i 21+4 



Pi 



~cm I 
Pi I Mi 



(P 



cm 



{p 



,cm-j2 



Xi 



1+2 



The virtual-photon width was defined by: 



r7 = r, 
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UK™\ Mt ) 2 +X 2 


K cm \ Mi _ 




(K cm ) 2 + X 2 



(16) 



(17) 



4 



A. Constraints on the Q 2 dependence for ol 

There were several physics constraints on the Q 2 dependence of the cross section which were realized in the fit form. 
For instance, the longitudinal cross section must vanish at the photoproduction point due to current conservation, and 
at large Q 2 due to helicity conservation in scattering from spin- 1/2 fermions. These constraints were independently 
imposed in both the resonant and non-resonant contributions to the cross section. This was accomplished for the 
resonant contribution by a suitable form for the longitudinal transition amplitudes, which vanished at Q 2 — » and 
Q 2 -> oo. 



B. Resonances included in the fit 



Seven possible resonance contributions were included in the fit to <tt, and represent those with the largest photo- 
couplings to the proton as listed by the particle data group [10| . These consisted of the first prominent resonance, 
the A P33(1232), two states in the second resonance region, the 5ii(1535) and the Z?i3(1520), and two in the third 
resonance region, the 515(1650) and the F 15 (1680). In addition, the Roper, J > 11 (1440) and a broad resonance around 
W ~ 1.9 GeV were included. Only the pion, eta, and 2-pion decay modes were included for each resonance and 
the branching fractions used are given in Table HI The region at W « 1.9 GeV is occupied by many states and 
was assigned an angular momentum of I = 3 in the fit. We have not tried to be exhaustive by including all known 
resonant states, but have, rather, tried to include enough of the dominant states to provide for a good fit to the 
data without requiring a burdensome number of additional parameters. For recent reviews of electromagnetic meson 
production in the nucleon resonance region we refer the reader to [Tl| . and to the latest fit to the unitary isobar 
model, MAID2007 [lj|. 

The resonances included for <tl were largely the same as for cry. However, due to the decreased sensitivity of the 
cross section data to <jl it was found that including more than a single resonance in the region of the -Di3(1520) and 
Sii(1535) around W ~ 1.5 had little impact on the goodness of the fit. Therefore, the second resonance listed in 
Table|I] (5x1(1535)) was not included in the fit for ox. The choice of excluding this particular resonance was arbitrary 
as it was found that excluding the Di3(1520) instead had a neglible impact. 





State 


ft* 


ft* 


ft 


1 


^33(1232) 


1.0 


0.0 


0.0 


2 


5u(1535) 


0.45 


0.10 


0.45 


3 


Di 3 (1520) 


0.65 


0.35 


0.0 


4 


Pis (1680) 


0.65 


0.35 


0.0 


5 


515(1650) 


0.4 


0.5 


0.1 


6 


Pn(1440) 


0.65 


0.35 


0.0 


7 


(1 = 3 assumed) 0.5 


0.5 


0.0 



TABLE I: Resonance number /, name (and quantum numbers), and branching fractions for the resonant states included in the 
fit. 



C. Resonance transition amplitudes 



For the transverse resonance transition amplitudes the fit form utilized was 



^' Q2 >= (i4voW^ 1 + TfwJ' (18) 



while for the longitudinal resonance transition amplitudes the fit form utilized was 



For large Q , A l T (Q ) reduces to 



^<« 2 > = (TT#$5IF-t' < 20) 
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which is just the dipole form for Cj = 2. For Q 2 — > 0, A l T (Q 2 ) linearly approaches the photoproduction value, Aj,(0), 
while A l L (Q 2 ) linearly vanishes. 



D. Non- resonant background 



For the non-resonant part of the transverse cross section the fit form which was found to describe the data well was 

\ - (°) 

^ (Q2 + a T^[bJ+cjQ 2 +djQ 



2 NR,i (n x 

«¥" = ^y. .„ ^3L^ n ,A Awy*''\ (2i) 



where 



with to w the pion mass, Qq = 0.05 GeV 2 , and AW — W~ m^. Previous resonance region fits have also included 
terms to describe the W dependence with (AW) l+1 ^ 2 . However, these did not attempt to fit beyond W ~ 2 GeV, 
where, in contrast to the data, this particular form continues to rise with increasing W 2 . Here it was found that this 
rise can be moderated by the additional x' factor. This allows a good fit to the data to be obtained up to W = 3.2 
GeV. The factor x' w x — Q 2 /2M p v for W 2 — M 2 » Mp-m^ and Q 2 » Q§, but has the important properties, 
required to fit the entire resonance data set, that it approaches a constant for Q 2 — > and vanishes at pion threshold. 

For the longitudinal cross section it was found that the non-resonant background could not be well described by 
powers of W 1 ^ 2 at fixed Q 2 , but was better represented away from pion threshold by (1 — x) a . Therefore, the form 
utilized was 

where 

x log(log(^ 



is a slowly varying function of Q 2 . For <r[^ R , we used Ql = 0.125 GeV 2 and m = 4.2802 GeV. 



(24) 



IV. THE DATA SETS 



A summary table of the data sets, including the number of data points and the Q 2 range of each set, is provided 
in Table HH 

The nucleon resonance region (W < 2 GeV) ep scattering data included in the fit came principally from a series 
of Jefferson Lab Hall C experiments which utilized the well- understood High Momentum Spectrometer (HMS). The 
various Hall C data sets have been estimated to have normalization uncertanties of w 2% and were found to agree 
in their regions of kinematic overlap to better than 1%. There also exists a large body of inclusive ep data from the 
Hall B CLAS collaboration in the kinematic region fit. However, in additon to having significantly larger systematic 
uncertainties, only results for the F% structure function have been published. Therefore, this data was not utilized 
since only cross section data were fit. 

For 0.18 < Q 2 < 4.5 the inclusive proton data utilized was dominated by the recent E94-110 high precision data 0, 
which have typical point-to-point uncertainties in e of less than 2%. This experiment was optimized to separate gl and 
ot from Rosenbluth separations [13] and therefore, provided a significant e range at fixed W 2 and Q 2 for constraining 
the separated cross sections. For 4.5 < Q 2 < 7.5 the resonance data aremainly from experiment E00-116 [14]. For 
this data set the uncertainties were typically statistically dominated. In addition, the cross sections at each W 2 and 
Q 2 were only measured at a single kinematic setting and, hence, these data did not provide any e range to constrain 
ol and <tt when taken by itself. Preliminary data from E00-002 [l5[ were used to extend to Q 2 range of precision 
data down to 0.05 GeV 2 . Data in the DIS region [l6[ (W > 2) were used to constrain the fit up to W = 3.1 GeV. 
Several sets of photoproduction data were used to constrain the fit at Q 2 = GeV 2 : data taken prior to 1975 [l7l - fl9T | 
and the more recent data from DAPHNE [2(|. While the latter is much more precise, it does not extend beyond the 
second resonance region. 
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Data Set 




Q Max 


# Data Points 




(GeV 2 ) 


(GeV 2 ) 




E94-110 [5] 


0.18 


5 


1259 


-bjUU-llo 1 14 1 


3.6 


7.5 


256 


E00-002 [15] 


0.06 


2.1 


1346 


SLAC DIS [16] 


0.6 


9.5 


296 


Photoproduction (Old) [17-49] 








242 


Photoproduction (DAPHNE) [20] 








57 



TABLE II: Summary of data sets included in the fit. 



CO 



g 




3000 - 



2000 - 



1000 



80 

60 

40 

20 


1.5 



E = 2.238, 



31 .93 




3.1 18, 9 = 61 .95 



0.5 - 




W 2 [GeV 2 ] 



W 2 [GeV 2 ] 



FIG. 1: (Color online) Comparison of the fit results (solid lines) to the JLab E94-110 resonance region data 5, 6] (solid triangles) 
versus W 2 , for several representative kinematic settings. Beam energies E are in GeV, and electron scattering angles 6 are in 
degrees. 



V. RESULTS AND DISCUSSION 



The results for the fit parameters are given in Tables ITU and ITVl Overall, the fit presented here provides an excellent 
description of the resonant structures seen in inclusive ep cross sections, as illustrated in Fig. [TJ for representative 
kinematic settings of the E94-110 experiment [f|. The fit does a reasonable job in describing the photoproduction 
data [l7l - [2Cjj (see Fig. [5]) , although it was not possible to obtain a perfect description of the dip region between the 
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W 2 [GeV 2 ] 



FIG. 2: (Color online) Comparison of fit (solid line) to the photoproduction data used in the fit (see Table Hi]) . 



I 


Mi 


r< 


A T (0) 


a» 


bi 


a 


Ai(0) 


di 


e. 


1 


1.230 


0.136 


7.780 


4.229 


1.260 


2.124 


29.4140 


19.910 


0.226 


2 


1.530 


0.220 


6.335 


6823.200 


33521.000 


2.569 


0.0 






3 


1.506 


0.083 


0.603 


21.240 


0.056 


2.489 


157.9200 


97.046 


0.310 


4 


1.698 


0.096 


2.330 


-0.288 


0.186 


0.064 


4.2160 


0.038 


1.218 


5 


1.665 


0.109 


1.979 


-0.562 


0.390 


0.549 


13.7640 


0.314 


3.000 


6 


1.433 


0.379 


0.0225 


462.130 


0.192 


1.914 


5.5124 


0.054 


1.309 


7 


1.934 


0.380 


3.419 


0.000 


0.000 


1.000 


11.0000 


1.895 


0.514 



TABLE III: Fit parameters for each resonance /, as defined in the text. Units of cross section are fib and all masses, momenta, 
and energies are in units of GeV. In addition to these parameters, the delta damping factor Xi was determined from the fit. 



A(1232) and the second resonance region. This is illustrated more clearly in Fig. [3l where two significant oscillations 
around unity are observed at low W and Q 2 < 0.1 GeV 2 (dominated by the photoproduction data). This is likely 
due to the Q 2 dependence for the transverse transition form factors chosen. It is certainly true that the individual 
transverse transition form factors in the second resonance region are not consistent with those extracted from exclusive 
analysis [11], although there is consistency in the overall transverse resonance strength in this region. Part of this 
inconsistency could be caused by the chosen form for the Roper amplitude A\. Recent data [Tl| indicate that this 



L/T 


i 


a L/T (°) 


a. 


h 


a 


di a 


T 


1 


246.1 


0.0675 


1.3501 


0.1205 


-0.0038 


T 


2 


-89.4 


0.2098 


1.5715 


0.0907 


0.0104 


L 


1 


86.7 


0.0000 


4.0294 


3.1285 


0.3340 4.9623 



TABLE IV: Values for the transverse (T) and longitudinal (L) non-resonant fit parameters. The units of cross section are 
fib and all masses, momenta, and energies are in units of GeV. In addition to those listed, the parameters Qo and m which 
determine the transition of a^ R to Q 2 = were determined from the fit. 
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amplitude decreases to zero by Q 2 s=s 0.4 GeV 2 and then increases again with increasing Q 2 . Such a Q 2 dependence is 
not possible with our chosen form for Aj, , and this could force a shifting of strength between the various resonances 
in this region which is Q 2 dependent. This is a topic that will be pursued in future fitting studies. 

The quality of the fit is quite good, overall, with nearly all data points differing from the fit by less than 5%, and 
more than half of the data points deviating from the fit by less than 3% as evidenced by Figure [5] As illustrated 
in Fig. 21 the fit describes data both at low e and high e reasonably well. This is especially true for Q 2 < 2 where 
the available data show little or no variation in the goodness of the fit for differing e values, indicating that both the 
transverse and longitudinal cross sections are well represented here. The data from E94-110 clearly indicate that the 
resonance structure in the longitudinal cross section becomes enhanced relative to the background for Q 2 between 
2-3 GeV, and that the resonanance peaks prefer widths that are narrower than for the transverse cross section. This 
indicates that the description of the data could be improved by allowing the resonance width parameters to be varied 
independently for the longitudinal cross section. 

The largest deviation with e is in the dip region between the second and third resonance region for 2 < Q 2 < 4 
GeV 2 . This is mainly due to the difficulty in fitting the W 2 dependence of the longitudinal cross section here. This 
region is dominated by the E94-110 data |5J which indicate significant resonance structure in the longitudinal channel 
in this Q 2 range, but with the longitudinal strength dipping very close to zero in the dip region [5] . This is inconsistent 
with a fit form for the non-resonant cross section which decreases monotonically as pion threshold is approached, and 
could indicate that our assumption of incoherency between the resonant and non-resonant scattering is starting to fail. 
Due to the lack of low e data at high Q 2 and W < 2 GeV, the fit for ax does rely to some extent on the extrapolation 
in Q 2 of the fit form used. 

In summary, we have developed a fit to the total inclusive transverse and longitudinal proton cross section that 
describes existing data at the 3% or better level over almost the entire range < Q 2 < 8 GeV 2 and 1.1 < W < 3.1 GeV, 
corresponding the to the kinematic settings available at Jefferson Lab with a 6 GeV beam. The fit can therefore be 
used to reliably evaluate radiative corrections and to extract spin structure functions from asymmetry measurements. 
The fit also provides a convenient representation of world data that can be used for the experimental evaluation of 
the high-x contribution to sum rules involving integrals over proton structure functions. 

FORTRAN computer code embodying the fit described in this article is available by email request from the authors. 
Tables of cross section data fit can be found in the relevant references, except for the preliminary data from JLab 
E00-002 which will be made available from the Hall C website at www.jlab.org/resdata. This webpage will be utilized 
as a repository for the final cross sections, as well as all available resonance region cross section data and fits. 

We thank V. Tvaskis for compiling the preliminary data table for Ref. [ll|. This work was supported in part by 
research grants 0099540 and 9633750 from the National Science Foundation. The Southeastern Universities Research 
Association (SURA) operated the Thomas Jefferson National Accelerator Facility for the United States Department 
of Energy under contract DE-AC05-84ER40150. 
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FIG. 3: (Color online) Ratio of all fitted cross section data to the fit results at the eight Q 2 ranges indicated. The units of 
Q 2 are GeV 2 . 
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FIG. 4: (Color online) Ratio of selected cross section data to the fit at the seven Q 2 ranges indicated (Q 2 in units of GeV 2 ). 
The data for e < 0.5 are shown as the solid triangles, while the data with e > 0.75 are shown as the open circles. 
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FIG. 5: Frequency distribution for the percent differences of the data to the fit for Q 2 > 0. 



